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The reported spatial variation in the fine-structure constant at high redshift, if physical, could 
be due to the presence of dilatonic domains, and one or more domain walls inside our horizon. An 
absorption spectrum of an object in a different domain from our own would be characterized by a 
different value of a. We show that while a single wall solution is statically comparable to a dipole fit, 
and is a big improvement over a weighted mean (despite adding 3 parameters), a two- wall solution 
is a far better fit (despite adding 3 parameters over the single wall solution). We derive a simple 
model accounting for the two-domain wall solution. The goodness of these fits is however dependent 
on the extra random error which was argued to account for the large scatter in most of the data. 
When this error is omitted, all the above solutions are poor fits to the data. When included, the 
solutions that exhibit a spatial dependence agree with the data much more significantly than the 
Standard Model; however, the Standard Model itself is not a terrible fit to the data, having a p- value 
of ~ 20%. 

PACS numbers: 



I. INTRODUCTION 

The universality of fundamental constants is one of the 
underlying tenets of physics. Of course some of these con- 
stants, like the fine-structure constant, a, may be dynam- 
ical and their values may be the result of an expectation 
value of some scalar field. Thus testing this universality 
is a bridge to physics beyond the Standard Model [![. 

Claims of a temporal variation in a from observations 
of quasar absorption spectra have been extended to in- 
clude a possible spatial variation as well @- The ini- 
tial indications of a temporal variation in a made use of 
the many-multiplet method Q and sparked an enormous 
amount of theoretical activity in attempts to explain 
it 0-11] • The Keck/Hires data which yielded a statisti- 
cally significant trend with Aa/a = (— 0.54±0.12) x 1CT 5 
over a redshift range 0.5 < z < 3.0 (the minus sign in- 
dicates a smaller value of a in the past). Subsequent 
studies based on VLT data using the same method have 
shown Act to be consistent with zero [t| [l(| • Of course 
these results can be made compatible if there is a spatial 
variation in a. 

Variations of a observed in individual absorbers using 
the many-multiplet method or more generally temporal 
variations, could be due to systematics related to the as- 
trophysical assumptions made regarding each absorber. 
For example, it is generally assumed that the abundance 
ratio of Mg isotopes take their terrestrial values. How- 
ever, even a slight enhancement in 25 ' 26 Mg/ 24 Mg could 
nullify many of the observed variations. Furthermore, 
such an enhancement in the heavier Mg isotopes could be 



explained by an earlier population of intermediate mass 
stars fill ] . 

Results from a recent VLT survey of 153 absorbers was 
performed 0], and taken together with the Keck data, 
leads to a large sample of 293 absorbers across the sky. 
However, it was shown, that there is a statistically sig- 
nificant spatial variation as demonstrated by a dipole fit 
to the data 



Aa 



m + A cos (i?) 



(1) 



which can not be accounted for by the systematic uncer- 
tainty due to the Mg isotopic abundances. In Eq. (JIJ, m 
is the monopole term, A is the magnitude of the dipole, 
and $ is the angle between the directions of the absorbers 
and the dipole. The template therefore has four parame- 
ters, m, A, and the two angles which specify the direction 
of the dipole. There is a marked improvement to the fit 
(using the full data set) using the dipole rather than a 
simple weighted mean (monopole) despite adding three 
parameters. 

If the fine-structure constant does depend on some dy- 
namical scalar field, cj>, it may well be space-time de- 
pendent. On cosmological scales, it is usually thought 
that the time variation dominates over spatial fluctua- 
tions, as suggested by most models. For example, sup- 
pose we couple a scalar field to electromagnetism through 
Bf(4>) p^pw ^ p^ v being the Faraday tensor and B F an 
arbitrary function of <f>. This will necessarily induce a 
coupling to matter which is generated radiatively if not 
present at the tree level (see below). The equation of 
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motion for the scalar field simply takes the form 
□0+^=0, 



(2) 



where V e g includes the self interactions of <j> as wen as 
any couplings to matter. If the Lagrangian contains a 
term Bn^ttinNN, then the coupling to matter is ef- 
fectively density dependent, and could serve as the source 
of spatial variations through 



n^ + m^tj) = B' N ((j))p N , 



(3) 



where is the scalar mass and pn is the baryon en ergy 
density, as is the case for the chameleon mechanism [12j . 
However, the density dependent shifts from the homoge- 
neous solution are typically extremely small except per- 
haps in the vicinity of a neutron star [1 31 - 41 5j j - In contrast, 
temporal variations are relatively easy to achieve partic- 
ularly over cosmological time scales, as long as the field 
remains light. 

Thus, we know of no physical or field theoretic model 
which could produce a dipole accounting for the observed 
spatial variations. Instead, we proposed [lj| to invoke the 
existence of a spatial discontinuity of the fine structure 
constant due to the existence of a domain wall crossing 
our Hubble volume (this idea was further studied in [l7|). 
In this case, a would take two values, a+ (the larger 
value of a) on our side of the wall and the second a_ at 
high redshifts on the other side of the wall. This implies 
that local constraints [l[ on the variation of a such as 
atomic clocks, Oklo and meteoritic dating will be trivially 
satisfied. 

The simplest way to implement this idea is to consider 
the following theory 



5' 



\MlR- - V{cjy) - \b f {<)>)F%, 



-gd 4 x, (4) 



where M~ 2 = 8irG is the reduced Planck mass. The 
scalar field (f> is assumed to have a simple quartic poten- 
tial 



(5) 



and a coupling to the Faraday tensor as well as to the 
fermions tpj. The coupling functions Bi are assumed to 
be of the form 



W)=-P^iJ- 1 + ^ (6) 

where the coefficients £j are constant and is a 
mass scale. This model depends on the parameters 
(A, M*, ?7, £p, £j) and we shall assume here that, at tree- 
level, only £f is non- vanishing. Nevertheless, the scalar 



field inevitably couples to nucleons radiatively through 
^ = m^(N\(^ F /4)F^\N) 0K This yields £ p = 
-0.0007£f and £ n = 0.00015£f 18] for the proton and 
neutron, respectively. Since most baryons in the universe 
are protons, we shall take £/v = £p for simplicity in our 
estimates. 

The main difference between the model studied here 
and previous models is that the scalar field is assumed to 
be heavy so that it is stabilized, hence we do not expect 
any local violation of the equivalence principle. Indeed, 
the current model does not exhibit any temporal varia- 
tion of constants once the phase transition has occurred. 
The resulting shift in the value of a on the other side of 
the wall is easily determined [l6[ 



a " M* 



(7) 



For simplicity, we shall assume r\ = M*, so that £p ~ 
1CP 6 is required to match the data. Then for A ^ 1 
and r\ ~ 1 MeV, one can show that such a wall makes 
only a tiny contribution to the overall energy density ( 



wall 



( 



•I 



) ) [3' ^ i s cosmologicallystable, and 



is compatible with microwave background 2C| and other 
astrophysical constraints. Furthermore, due to the dy- 
namics of the phase transition producing the wall, we 
expect to be left with order one large wall per Hubble 
radius [2ll | which is moving slowly towards us jl6] . 

Given the large data set available 0, y] , it is possible 
to be more quantitative concerning the wall. We can in 
fact use the data to determine the position of the wall, 
and the potential drop across it. We can further make 
a statistical judgement as to whether the wall is an im- 
provement over a simple temporal variation (a monopole 
fit), and compare it directly to the dipole fit found in 0. 
As we will see below, while in ideal circumstances, the 
wall is an approximation to a dipole, the direction of the 
wall only lines up well with the dipole when using the 
Keck data alone. For the full data set, best fit position 
of the wall is not aligned with the dipole. Nevertheless, 
the fit for the wall as determined by a \ 2 analysis is com- 
parable to that of the dipole. While less significant than 
the Keck data, the VLT data shows a tendency for pos- 
itive variations of a. We show that a two- wall solution 
(despite requiring three additional parameters) is a far 
superior fit to the data than either the dipole or single 
wall solution. 

In the next Section, we present a model with three 
dilatonic solutions, which is the basis for the two- wall 
fits of the data discussed below. This model is a simple 
extension of the model of [llj]. In Section Hill we briefly 
describe the data that we use. In Section HVl we outline 
the key algorithm used in the walls fit, used to deter- 
mine in which vacuum an absorber lies. In Section [V] 
we present the results of our fits. The significance for 
a spatial dependence in the data, and the orientation of 
the dipole and wall fits are discussed in Section I VII Our 
conclusions are summarized in section fVIII 
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II. THE TWO- WALL MODEL 



In many extensions of the Standard Model, the value 
of a is a function of the dilaton field, as indicated in (M, 
In general, <f> is the real part of a complex field In [la ], 
it was suggested that the potential of <f> has two minima, 
and that we live in one of them, while some of the ob- 
served absorbers live in the other vacuum. In this paper, 
we also fit the data against a two- wall model. In this fit, 
absorbers have either the same value ckq that we mea- 
sured on Earth, or the value ao (1 + A), or ao (1 — A). 
In this Section, we briefly indicate how the model of [la] 
can be modified to have three dilatonic vacua. 

The relevant term in the action are 



C = -\d^-V{$)-\B F {<t>)Fl v 



V 



a m 



2 



b = Re$ 
V , (8) 



where the three parameters A, 77, e are positive and real 
numbers. We note that A is dimensionless, while 77 and 
e have mass dimension one. The constant Vq is chosen 
such that the potential vanishes at its minimum. For 
computational simplicity, we assume that e <C A 77, so 
that the second term in the potential can be treated as 
a perturbation. The potential has the three minima 
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V 2 



n = 1, 2, 3 



9e 2 3e _ f j_ 



(9) 



The phases of the three minima correspond to the com- 
plex vectors shown in the right half of Figure [TJ We as- 
sume that we live in the n = 1 vacuum. To evaluate the 
values of a in the other vacua, we assume that B is still 
given by eq. ©, with 4> — Re and expand it for small 
£. We obtain 
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(10) 



For e <C A77, the radial excitations of <f> around the min- 
ima are much more massive than the angular excitations. 
At energies below the mass of the massive radial excita- 
tions, we have a consistent physical description which in- 
cludes only the light degrees of freedom. At leading order 
in e, the canonically normalized field, 8 C , corresponding 
to the light excitations is 



d> = -?Le^+0(e) 
and is controlled by the lagrangian 
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C = --(d6 c ) 2 -erf 
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FIG. 1: The left side of the Figure illustrates a cosmic string 
formed at the first stage of symmetry breaking; the arrows 
indicate the complex value taken by $ in different spatial po- 
sitions. The central dot represents the string (which extends 
out of the page). At the second stage of symmetry breaking, 
<E> settles in the different minima inside separated domains 
as seen in the right side of the Figure. The solid lines repre- 
sent sections of the walls separating the domains. This part of 
the Figure also illustrates a typical wall configuration in our 
fits: our location is at the center of the dotted circle, which 
represents the typical radial distance to the Keck and VLT 
absorbers. 



In the following, we assume A — O (1) for definiteness. 
At high temperatures, the potential © is modified by 
thermal effects. At temperatures greater than O (77) the 
model is in the unbroken phase, namely the potential is 
minimized at $ = 0. Below this temperature, the nearly 
U(l) invariant set of minima |$| ~ appears. This 

residual U(l) symmetry is preserved as long as thermal 
fluctuations are energetic enough to cross the potential 
barrier between the three vacua ([9]); this is the case until 

T ~ O ((jj^j v\- Below this temperature, the system 

is in the fully broken phase. 

Therefore, for e <C 77, symmetry breaking occurs in 
two stages. Cosmic strings are formed at the first stage. 
Outside each string, $ has fixed magnitude and variable 
phase. As the temperature decreases, domains of the 
three vacua ([9|) form; the different domains are separated 
by a domain configuration with |$| = + O (e), and 
with a phase interpolating between the values in the two 
domains. As illustrated in Figure [TJ at least three domain 
walls stream off each string (more domains will be present 
for strings characterized by a higher winding number). 

For a planar and static wall perpendicular to the z 
axis, it is convenient to define 
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T 



erjz, 



(13) 



so that the action is manifestly put in the form of the 
sine-Gordon model. Indeed, starting from (|12[) . and per- 
forming these redefinitions, we obtain the domain wall 
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action 
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dtdxdy 



cos 6 + 1 



dz 




V 



(14) 



In these coordinates, the three minima © correspond to 
the 9 — — 7r, 7r, 37r minima of V, respectively. The wall 
solution is obtained from a standard BPS procedure. The 
equation of motion ^ — = is integrated to give 

£*h£ (lu 




V 



0. 



(15) 



This equation can be further integrated to give (returning 
to the original coordinates) 



'/ 



tan 



(16) 



This is the domain wall solution interpolating between 
the two minima n = 1 (at z = — oo) and n — 2 (at 
z = +oo) of ©. We note that the thickness of the wall 

is of 

The tension of the wall is most easily computed using 
the rescaling (fT3|) . and using (TTBl . we obtain 



T = 



energy 



= dz 




(17) 



The cosmological and astrophysicaHimits for the single 
wall model were discussed in Ref. [l6j]. This discussion 
can be readily extended to the current case. For example, 
the requirement that the CMB is not distorted by the 
walls translates into an upper limit [161 ] 



1/6 



■q < few MeV , 



(18) 



on the wall tension. The tension of the strings is also of 
order ij and values compatible with (jT8"f do not lead to 
any CMB limit from the strings. 

Another relevant constraint is related to avoiding ex- 
cessive emission of dilatonic quanta in supernovae. In the 
current model, such quanta in the vacuum n — 1 have a 
squared mass ml — 9 e r\ (these are the light excitations 
described by ([T2|) ). If this mass is smaller than the su- 
pernovae temperature, a large number of dilatonic quanta 
are produced in supernovae. Such quanta decay back into 
photons with a rate T ~ ^rnh/Ml ~ £| e 3 / 2 if/ 2 /Ml . 
If the corresponding decay length is smaller than the size 
of the supernova core, we do not have energy loss into 
such quanta [l6| . This occurs for 
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e > O (1(T 2 ) MeV 

1 & J V 



(19) 



(for typical supernovae temperatures of order T ~ 
30 MeV). If we make the natural choice r\ — O (M*), we 
see that the two conditions (|18|) and (fT9|) indeed allow 
the hierarchical choice e <§; fy. 

Next, we need to discuss the modification of the poten- 
tial due to the interaction of $ with matter. As discussed 
in [16| , the interaction generates an additional term in the 
potential 



at, ^ Re$ 

Ay = 6V -jrjr- Pb ■ 



(20) 



where is the energy density of baryons in the universe, 
and £at is the constant controlling the coupling of the 
dilaton to nucleons in ([6]). As in [lj|, we assume that 
this coupling is generated radiatively starting from the 
coupling of the dilaton to photons [7j, and we disregard 
the contribution of neutrons. This gives £at — — 0.0007 
[HI]. We then find, for the three vacua ([9]) 



AVi = , AV 2 



4 x 10" 



= +A fc 

io- 



, AV 3 = - 



(21) 



where p^ 



-48 



GeV is the current value of the 



1.8x10" 

baryon energy density. 

The coupling to baryons breaks the degeneracy be- 
tween the energies of any two domains separated by a 
wall. This generates motion of a wall towards the re- 
gion of greater potential energy, however, the motion is 
slowed down by Hubble friction. The velocity v of the 
wall is described by [l6| 



at tension W aii 



(22) 



where R is the radius of the universe and 7 = 1/ yl — v 2 . 
For the single wall model studied in [l6[ , this resulted in 
a current wall velocity, vo ~ 0.004. The only significant 
difference between the single and the double wall model is 
the parametric dependence on e <g. rj of the wall tension. 
This results in an increase of O (y/f) for the value of 
in the model considered here with respect to the one in 
16]. This means that for the model with two walls, we 
estimate vo ~ 0.004 as long as the current motion is 
non-relativistic. The two bounds (|Tg)) and (|T9")) suggest 
a hierarchy e/rj ~ 10~ 2 , in which case it is also safe to 
disregard the motion for the two-wall models. 

A second potential worry associated with (|2~T1) is the 
tunneling of the domains n = 1, 2 into the domain of 
lowest energy, n — 3. The tunneling rate is suppressed 
by [27|,[2| 



r oc exp 



27tt 2 T 4 



2 Ay 3 



exp 



10 



13") 



(1+*)- 



(23) 



where T is the tension evaluated in (ITT)) and AV is the 
potential difference (f2Tj) . In the numerical evaluation, we 
have set e = 10~ 2 ?/, rj = M» = MeV, £ F = 10~ 6 . The 
wall forms at z ~ 10 9 . We see that the tunneling rate is 
always negligible. 
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III. DATA USED 

The available data on possible variations in a come 
from two sources: Keck data 3] containing 140 absorp- 
tion systems and VLT data Q containing 153 systems. 

Concerning the VLT measurements, we use the re- 
sults appearing in Table Al of the second reference in 
2]. Consistently with what done in 0], we discard the 
fourth absorber in the list as an outlier. We use the 
last column of that table for the values of Aa/a for 
each of the absorbers. In several fits using these data, 
ref. Q added in quadrature a common "random error" 
o'ran = 0.905 x 10 -5 to the value of a given in the table 
for each absorber. The motivation given in [2j is that 
the VLT data are too scattered, and therefore the origi- 
nal error given in the table must be underestimating the 
true error. Ref. |2| actually uses slightly different values 
for er ran for fits of the data to different functional forms. 
The mathematical procedure leading to the values used 
is described in their Section 3.5.3. This procedure ulti- 
mately allows the data to agree with the fitting model 
(the random error is determined by requiring that a sub- 
set of the data - once the biggest outliers are excluded - 
has a reduced x 2 equal to one for the model under which 
the data are being fit). This is certainly a necessary con- 
dition for a meaningful model comparison, if one wants 
to claim that the better model fits the data. However, we 
stress that the procedure is arbitrary to a large extent, 
and that ref. Q did not identify the physical source of 
this additional error. We therefore regard as unmotivated 
any p-value obtained for the fits when cr ran is added; we 
agree that it is nonetheless interesting to compare the 
significance of different fits to the data, in the hope that 
the physical source of the error is eventually found, and 
that this error is not correlated with the spatial position 
of the observers. Concerning this last comment, we stress 
that, for er ran = 0, all of the models considered in [2j and 
in this paper are a very bad fit to the data, and therefore 
we cannot rule out that any claimed spatial dependence 
is spurious. In our data analysis, we compute the p-value 
of the fits with and without adding £r ran , in order to see 
how the significance of the dipole and wall fits is affected 
by this random error. 

Concerning the Keck/Hires data, we use the same data 
studied in . Such data were analyzed in [H, [24| , and 
can be found in the On-line Data Table 1 linked to [25| . 
This table contains 143 absorbers, 27 of which are de- 
noted as a "high-contrast sample"; such absorbers are 
at redshift z > 1.8 and both strong and weak transition 
lines had to be fit to extract the value of a associated to 
them. According to [23[ this is the origin of the excess 
scatter present in these data, which is not adequately re- 
produced by the calculated statistical error. To account 
for this, ref. 23[ added in quadrature the common error 
<T = 1.75 x 10~ 5 to the individual error of each high- 
contrast absorber. The value of this common error was 
chosen in [23[ so that, once the total error is used, a com- 
mon mean fits these 27 absorbers with a reduced x 2 equal 



to one. Ref. Q used different values for this common er- 
ror according to different fits of the data, analogously to 
what done for the VLT data. In our study, we compare 
the significance of different models both with and with- 
out the additional error, a — 1.63 x 10 -5 , which is the 
value used in 0] for most of the dipole fits. We stress 
that this error is added only to the 27 "high-contrast" 
absorbers. According to ref. [26|, the original values of 
a studied in [23| are affected by an incorrect wavelength 
calibration. Ref. [24j estimated the shift induced by 
the miscalibration on each individual absorber. In most 
cases, the quoted values are smaller than the statistical 
error on each measurement, and ref. [24| concluded that 
this effect does not significantly impact the data analysis. 
This recalibration was disregarded in the analysis of Q ■ 
Additionally, ref. |2J indicated that two absorbers in the 
table in [25j, and included in the previous analysis, are 
affected by overly large calibration problems, and should 
be disregarded. One additional absorber is also identified 
as an outlier in the analysis of @, and disregarded. To 
compare our results with those of 0, we also disregard 
the recalibration and these three absorbers in our study. 



IV. WALL FITTING PROCEDURE 

Let us first discuss our fitting procedure for the sin- 
gle domain wall model. The wall separates the dilatonic 
vacuum we live in (characterized by the value of a mea- 
sured on Earth) from the other vacuum (characterized 
by a different value of a). We need to establish an algo- 
rithm that, for any wall configuration, indicates whether 
a given absorber is in our vacuum, or in the other one. 

As discussed above, the wall can be assumed to be 
planar and static today. This means that one can choose 
coordinates so that the wall is at a constant comoving 
cartesian coordinate, say x = x wa n . The wall has a point 
which is closest to us. Let us denote this point by P. 
The position of P is characterized by the two angular 
polar coordinates (#*, </>*) and by the redshift z*. An 
additional parameter of the model is the value of a in the 
other vacuum. These four parameters completely specify 
the wall model for the purpose of data fitting. 

The dashed line shown in Figure [5] illustrates how the 
wall appears on our sky. In the Figure, the coordinates 
have been chosen so that the origin corresponds to our 
location, while points at greater radial distance from the 
origin correspond to points at greater redshift on our sky 
(we did not attempt to show this relation to scale in the 
Figure). The coordinates in the Figure have also been 
chosen so that the point P lies on the positive vertical 
axis, = 0. The wall is described by a function 8 (z), 
that relates the redshift z of a generic point on the wall 
to the angle 8 formed by the line of sight to this point 
(as seen by us), and the line of sight to the point P. To 
obtain this function, we first note that the radial comov- 
ing coordinate of this generic point on the wall is related 
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wall 




integration variable z' 



1, this relation becomes 



FIG. 2: Illustration of how we determine whether an absorber 
is in our vacuum, or in a different one. Our location is placed 
at the origin, and increasing radial distance correspond to 
increasing redshift. The Figure shows the wall as it appears 
on our sky and an absorber at redshift za- This redshift is 
shown as the dotted circle. The angle 9 (za) is the angle at 
which we see the wall at that redshift; the angle 9wa is instead 
the angle between the center of the wall and the absorber. 
The case shown has 8wa > 9 (za) and the absorber is in our 
vacuum. 



to this angle by 



cos6> 



(24) 



where is the radial comoving coordinate of P. We 
then recall the relations adr = dt = da/ (aH), where a 
is the scale factor, and H the Hubble rate. We insert 
the expression of the scale factor in terms of redshift, 
a = ao/(l + z) in this differential relation (ao denotes 
the value of the scale factor today) and we integrate the 
resulting equation to obtain a relation between r and z. 
We finally insert this result into (|24l) . to obtain 
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(25) 



where we have restricted our attention to a flat uni- 
verse containing only matter and a cosmological con- 
stant. f2A,o — 0.728 and Q m ,o = 1 — ^A,o denote the 
present fractional energies of cosmological constant and 
matter, respectively (the numerical value is the Maximal 
Likelihood value given in [22j for a ACDM universe, us- 
ing WMAP7, BAO, and H data). Upon changing the 
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F{z*) 
T(z) 



dy 



y 



1— Oa,o 

(26) 

This relation correctly reproduces 9 (z* ) = for the 
case in which the generic point on the wall coincides 
with P. It also indicates that the wall is seen within 
a maximal opening angle, which is mathematically given 
by 9 (z = oo) (for fitting purposes, there is no difference 
between considering the maximal opening angle given by 
z = oo and the one given by the redshift at which the 
wall formed, since the latter quantity is much greater 
than the redshifts of the absorbers) . The opening angle 
is maximal for z* — 0, when the wall passes through our 
location. In this limit eq. (f2l)f gives 9 — n/2, which 
correctly indicates that the wall is seen in half of the sky. 

It is now straightforward to describe the algorithm that 
we use to determine whether an absorber is in our vac- 
uum or in the other one. Figure [2] also illustrates the po- 
sition of a generic absorber. For each wall configuration, 
and for each absorber, we compute the angle between 
the line of sight of the absorber and the line of sight of 
P. This is the angle 6>wa indicated in the Figure. We 
then insert the redshift of the absorber za into eq. (|26l) . 
The corresponding quantity 9 (za) is also shown in the 
Figure, and it corresponds to the angle between the line 
of sight of P and the line of sight of a point on the wall 
having the redshift za- If, as in the Figure, 6*wa > 9 (za), 
the absorber is seen by us outside the region delimited 
by the wall for that redshift. Therefore the absorber is 
in our vacuum. If instead 6*wa < 9 (za), the absorber is 
"beyond the wall" , and therefore in the other vacuum. 

In the case of two walls, one could have the situation in 
which (i) the walls are far apart, so that they do not in- 
tersect (at least in the region occupied by the absorbers), 
or (ii) they are close enough that they do intersect. In 
case (i), the stationary and planar approximation holds 
for both walls, and we readily extend to this case the pro- 
cedure just outlined for the single wall case. In the case 
(ii), one should in principle solve the field theory equa- 
tions at the intersection of the two walls to determine 
the precise spatial distribution of the different vacuum 
domains. For simplicity, we disregard the configurations 
of the type (ii) in our analysis. In practice, referring to 
Figure [TJ we imagine that the string connecting the walls 
is at a comparable but slightly greater distance relative 
to the Keck and VLT absorbers. Concretely, in our fits 
of the two-wall model, we first choose the 6 parameters 
that characterized the geometry of the walls (three per 
wall, as described above), implicitly assuming that the 
walls do not intersect. We then use (in our numerical 
fitting program) the algorithm described above to deter- 
mine whether an absorber is beyond either of the two 
walls. If the algorithm gives that no absorber is beyond 
both walls, then the assumption that the two walls do 
not intersect (in the region occupied by the absorbers) 
is correct, and the configuration is indeed of the type 
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(i). Otherwise, the configuration is of the type (ii), and 
we disregard the initial choice of the 6 parameters in the 
data analysis. Therefore, when we fit the two-wall model, 
we are only fitting the data against a subset of possible 
configurations. This effectively results in the fact that 
the significance that we quote for the two-wall model is 
a conservative figure: we cannot exclude that a configu- 
ration of the type (ii) would be a better fit to the data 
than those that we probe in our analysis. 



V. RESULTS 

We perform separate fits for the Keck measurements 
alone (140 objects), for the VLT measurements alone 
(153 objects), and for the total combined set (293 objets). 
The data were described in Section IIIII The data are 
given in Equatorial J2000 coordinates; the polar angles 9 
and <fi correspond to "declination" (given in degrees) and 
"ascension" (given in hours), respectively. The relation 
between the different units is 



9 7r / declination 

radians ~~ 2 V 90° 



7r ascension 



radians 12 hours 
The data are fit against the 



(27) 



• Standard Model: no variation of a with respect to 
that measured on Earth (0 free parameters); 

• a monopole template: a common constant m = 

, where a$ is the terrestrial value (1 free pa- 
rameter) 

• the dipole template ([T]) (4 free parameters) 

• the one wall model of [16[ (4 free parameters) 

• the two- wall model introduced in Section [TT] (7 free 
parameters) 

Two parameters of the dipole fit are the dipole and 
monopole amplitudes A and m specified in (p}. The 
other two parameters, 9 and 0, are the angles that specify 
the direction of the dipole in a given coordinate system 
(specifically, we use Equatorial J2000 coordinates in all 
our fits). The angle 9 should not be confused with the 
angle # defined after eq. ((TJ). The two angles 9 and <\> 
actually point towards a single position on the sphere, 
and therefore also specify the orientation of the dipole. 
Reversing this orientation (9 — s- tt — 9 , <fi — > <f> + ir, if the 
angles are expressed in radians), and changing the sign 
of A in ([1} results in the same dipole fit. 

The four parameters characterizing the one wall model 
are the two angular coordinates and the redshift of the 
point of the wall closest to us, and the value A = a ~°° 
corresponding to the vacuum on the other side of the 
wall with respect to us. In the two-wall model there are 



three additional free parameters characterizing the posi- 
tion of the second wall; in principle, an eighth parameter 
could be immediately introduced, since one can trivially 
modify the model introduced in Section [TT] to have two 
independent values Aj beyond the two walls. We remove 
half of the freedom by imposing that the values of A* 
beyond the two walls are opposite to each other. This 
is done to reduce the dimensionality of the parameter 
space, in order to facilitate the search for the best fitting 
parameters. Due to this choice, changing the sign of A 
and the order of the two walls results in the same config- 
uration. Finally, as mentioned at the end of the previous 
Section, we exclude from the fit the configurations for 
which the two walls intersect each other (at least, within 
the region occupied by the absorbers). This assumption 
is made for simplicity (otherwise we would need to study 
the dynamics of the intersection, or at least introduce 
additional parameters specifying the position of a third 
wall, see Figured]). 

For any fit, the x 2 of any parameter choice is computed 
in the standard way 



X 2 [parameters] 



(Act | _ Aa I \' 
\ a \i a I parameters/ 



(28) 



where the sum is performed over the absorbers in a given 
dataset, and where Ui is the error on the — |j mea- 
surement; — | parameters is instead the theoretical value 
for that parameter choice. The likelihood is defined as 

The dipole and walls fits are performed with a Markov 
chain. Each point in the chain corresponds to a point 
in the n— dimensional space of the n parameters that are 
being fit. The first point in the chain is chosen at ran- 
dom. The chain is then characterized by the algorithm 
that allows one to add the n + 1— th point to it starting 
from the n— th point. The addition is performed by iden- 
tifying a candidate point and by accepting it with a given 
probability. The identification is performed by doing a 
small random step in parameter space starting from the 
n— th point in the chain. The candidate point is accepted 
with the probability 



probability accepting point n + 1 = Min 



+ 1 



r 

*~"n 



(29) 

where C n is the likelihood of the n— th point, and analo- 
gously for n+1. We note that a candidate point is always 
accepted if it has a better likelihood than the point that 
was last accepted in the chain; otherwise, the probability 
decreases in proportion to how worse the likelihood of the 
candidate point is. A long chain then behaves as a grid 
which is more dense in regions with higher likelihood. 

The point in the chain with smallest % 2 is our solution 
for the parameters that best describe the data within a 
given fit. To find the confidence intervals around this 
best point, we vary the value of one parameter at a time, 
and marginalize over the other parameters (namely, the 
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Fit 


Keck 


VLT 


Combined 


Std. Model 


217.1 (3 x 10" 5 ) 


280.2 (2 x 10" 9 ) 


497.3 (10~ 12 ) 


Monopole 


180.3 (10~ 2 ) 


269.9 (10~ 8 ) 


496.5 (10 -12 ) 


Dipole 


176.5 (10~ 2 ) 


250.8 (4 x 10~ 7 ) 


449.6 (4 x 10" 9 ) 


One wall 


165.9 (4 x 10" 2 ) 


229.5 (3 x 10" 5 ) 


449.2 (5 x 10 -9 ) 


Two walls 


162.0 (4 x 10~ 2 ) 


209.0 (5 x 10~ 4 ) 


397.0 (10~ 5 ) 



TABLE I: x 2 (and p— values) without adding random errors 
on the VLT data or the 27 "high contrast" Keck data. 



other parameters are free to vary, until the configuration 
with the minimal x 2 is found). For the dipole fits, the x 2 
distribution is approximately Gaussian around the best 
fit point. Therefore, we determine the (approximate) la 
confidence level by varying each parameter until the x 2 
(with the other parameters marginalized over) increases 
by 1 with respect to the best fit point. For the wall fits, 
this was possible only for the parameter A and because 
the distribution does not resemble a Gaussian, the 68% 
CLs are determined directly from the the likelihood func- 
tion. For the other parameters, the x 2 behavior is too 
irregular. We show the x 2 distribution for some of the 
fits in Figures HH [3 and E 

The results of our analysis are summarized in a num- 
ber of Tables and Figures. In the Tables, we give the x 2 
and the corresponding p— values of a number of fits to the 

data. We recall that the p value p = T f§, *P) /T (§) 

(where v represents the number of degrees of freedom) 
indicates the probability that a set of fictitious data gen- 
erated from the model being studied has a greater x 2 
than the fit of the actual data. Therefore, it is a measure 
of how well the model fits the data. 

The results given in Table U are derived when omit- 
ting the additional random errors discussed in Section 
Hm For the Keck data, we see a significant drop in x 2 f° r 
the monopole which was the basis of the claim in [3j for 
a temporal variation in a. There is little motivation here 
for a dipole, but the single wall solution does bring an- 
other significant drop in x 2 ■ Here, adding a second wall 
does very little. For the one-wall solution, the p-value 
reaches 4 %. For the VLT data, the monopole solution 
offers a modest improvement in x 2 , and we see that x 2 
has significant improvements as we move from the dipole 
to the two-wall solution. In each case, the p-value is 
extremely small. For the combined data, the monopole 
offers almost no improvement, while the dipole and one- 
wall solutions give an almost identical drop in x 2 , which 
is further lowered in the two-wall solution. In total, x 2 is 
lowered by 100, at the cost of seven parameters. As we 
see from the very low p-values, none of the fits are good 
representations of the data. 



combined data, respectively, with the additional random 
uncertainty included. The basic behavior of x 2 described 
for Table U is repeated when the extra random errors are 
included. Of course, overall values of x 2 are significantly 
lower with correspondingly high values of p. We stress 
that the value of the extra error is determined a posteri- 
ori during the data analysis itself @, [23[ , and its value is 
chosen so that the dipole fit is a good fit to the data. 

We note that as seen in Fig. [3l there are multiple local 
minima for the value of A corresponding to the shift in 
a. For the Keck and VLT data alone, the lowest two 
minima have comparable values of x 2 and as such, there 
are two distinct ranges for A at the 68% CL: (-1.24 - 
-1.00) and (-0.97 - -0.68) for Keck, and (1.01 - 1.68) and 
(1.78 - 2.33) for VLT. For the combined data, the 68% 
CL falls within a single range for A. 



Keck Fits 


Parameters 


X 2 (p-value) 


Std. Model 




157.9 (0.14) 


Monopole 


m = -0.57 ±0.11 


132.8 (0.63) 


Dipole 


m = -0.47 ±0.15 
A = 0.41+g 

= -47t£, 0=16.Ot 2 ; 6 3 


131.0 (0.61) 


One wall 


A = -O.SOl^ 2 , z = 0.15 
6> = 32 , = 4.0 


127.6 (0.68) 


Two walls 


zi = 0.16 , 0i = 33.9 , 0i = 3.9 

22 = 0.74 , 02 = -37.6.5 , 02 = 18.2 

A = -0.82 


127.5 (0.62) 



TABLE II: Fits for the Keck data, adding random errors on 
27 "high contrast" data. The quantities m, A, A are given in 
units of 10" 5 . 



VLT Fits 


Parameters 


X 2 (p— value) 


Std. Model 




152.5 (0.50) 


Monopole 


m = 0.21 ±0.13 


149.8 (0.54) 


Dipole 


m = -0.11 ±0.19 
a _ i 17+0.47 

= -62 ± 14 , = 18.3+j-j 


141.8 (0.65) 


One wall 


A = 1.38l°;^ , z = 1.34 
= -37.8 , = 20.1 


132.7 (0.83) 


Two walls 


zi = 1.04 , 0i = -42.9 , 0i = 19.3 

22 = 1.50 , 02 = -14.6 , 02 = 8.1 

A = 1.44 


127.2 (0.87) 



TABLE III: Fits for the VLT data, including random errors. 



Tables imillll and II VI present the solutions for the best 
parameter values, the x' \ and the p— value of the vari- 
ous fits for the Keck data, the VLT data, and the total 



We see that the dipole or wall solutions offer little im- 
provement for the Keck data alone. The single wall gives 
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Comb. Fits 


Parameters 


X 2 (p— value) 


Std. Model 




310.4 (0.23) 


Monopole 


m = -0.22 ±0.08 


303.7 (0.31) 


Dipole 


m = -0.18 ±0.09 

A = 0.97 ±0.21 

6> = —61 ± 10 , (j> = 17.311;? 


280.6 (0.63) 


One wall 


A = -1.061^2 » z = °- 45 
6 = 68.3 , <\> = 20.2 


282.0 (0.60) 


Two walls 


Zl = 0.44 , 6»i = 67.9 , <f>i = 20.1 
z 2 = 1.02 , 6 2 = -39.4 , 4> 2 = 19.2 
A = -1.12 


263.2 (0.83) 



TABLE IV: Fits for the Keck and the VLT data, including 
random errors on 27 "high contrast" Keck data, and on all of 
the VLT data. 



310 
305 
300 
295 
290 



285 
280 




A * 10 



FIG. 3: x 2 vs - A for the one-wall model fit of the total data; 
see the main text for details. 



a large improvement for the VLT data. For the combined 
data, we see again that the single wall and dipole fits are 
clearly better than the Standard Model or the monopole 
fit, while the two-wall fit offers the best solution to the 
combined data. 

Figures [3j |H O El correspond to the one- wall fit of the 
combined data (with the additional random uncertainty 
included). In each Figure we show the x 2 obtained by 
keeping one parameter fixed (to the value shown on the 
horizontal axis of the Figure) and by varying the remain- 
ing 3 parameters until the minimum x 2 is found (i.e., by 
marginalizing over the remaining 3 parameters) . For each 
Figure, 200 values of the fixed parameter plus the result- 
ing x 2 are connected by a solid line. For each value, the 
marginalization procedure is done with a Markov chain 
of 500, 000 points. The large dots visible in three of the 
Figures are obtained with a chain of 20 million points, 
and verify that the shorter chain used to obtain the solid 
line is sufficient for the marginalization. 

In Figure[7l we show the position of the Keck absorbers 
on our sky and the wall configuration for the single wall 




FIG. 4: x vs - redshift of the closest point on the wall, for 
the one-wall model fit of the total data. 




-80 -60 -40 -20 20 40 
e [Degrees] 



60 80 



FIG. 5: x vs - declination of the closest point on the wall, for 
the one-wall model fit of the total data. 
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292 
290 
288 
286 
284 
282 
280 




10 15 
[Hours] 



20 



FIG. 6: x 2 vs - ascension of the closest point on the wall, for 
the one-wall model fit of the total data. 
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best fit given in Table |TTJ Each absorber is marked by 
a zero sign if it is in our vacuum (since Aa = in this 
case), of by a minus sign if it is in the other vacuum (since 
Aa < in this case). The location of the point on the 
wall which is closest to us is marked with an asterisk. The 
wall boundary, shown with a dashed line in the Figure, 
corresponds to the (mathematical) locations at which the 
wall would be seen at infinite redshift on our sky. The 
wall intersects any line of sight pointing in the region 
limited by that boundary (and including the asterisk); 
the redshift of the intersection is an increasing function 
of the angle between that line of sight and the line of sight 
of the nearest point of the wall (this corresponds to the 
angle denoted by 9 (za) in Figure[2|). An absorber in that 
portion of the sky can either be in our vacuum if it is "in 
front of" the wall (namely if its redshift is smaller than 
the redshift at which the wall intersects the line of sight 
of that absorber), or in the other vacuum if is is "behind" 
the wall. All the absorbers in the other portion of the sky 
(the one delimited by the dashed line boundary, and not 
containing the asterisk) are instead in our vacuum. Also 
shown in Figure [7] is the position of the dipole solution, 
marked by a 'D', which is well aligned (perpendicular) to 
the wall. 




FIG. 7: Keck absorbers (marked according to which vacuum 
they are in) and wall configuration for the single wall best fit 
given in Table [III The asterisk indicates the position of the 
point on the wall closest to us; the dashed line indicates the 
boundary of the wall. The letter D marks the direction of the 
dipole, oriented where the dipole amplitude is greatest. 

The VLT absorbers and the corresponding single wall 
best fit (Table UrT|) are shown in Figure [H The absorbers 
in the other vacuum are marked with a plus sign in this 
Figure, since Aa > in that vacuum. The wall solution 
in this case is found in a different region of the sky, though 
it too is well aligned with the dipole (which moved very 
little). The combined sample and the corresponding sin- 
gle wall best fit (Table llVf are shown instead in Figure |U 
In this case, the single wall attempts to interpolate be- 
tween the two previous solution and is no longer aligned 
with the dipole. Finally, Figure [10] shows how the com- 
bined sample separates into three vacua for the two- wall 
best fit (Table IIV|) . We note the presence on the sky of 




FIG. 8: VLT absorbers and wall configuration for the single 
wall best fit given in Table [TTT1 




FIG. 9: Total absorbers sample and wall configuration for the 
single wall best fit given in Table ITVl 




FIG. 10: Total absorbers sample and wall configuration for 
the double wall best fit given in Table ITVl 



a region where the two wall solutions intersect. The ac- 
tual wall solution should be modified in that region so 
to properly account for this (we expect the presence of a 
third wall, see Figure [T]). However, all absorbers in that 
portion of the sky have a redshift smaller than the inter- 
section (in fact, they have a redshift smaller than the one 
corresponding to the second wall). Therefore, the more 
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Statistical preference of 
Standard Model over : 


Keck 


VLT 


Total 


Monopole 


10" 6 


0.10 


0.012 


Dipole 


4 x 10" 5 


0.03 


7 x 10" 6 


One wall 


7 x 1CT 6 


3 x 10~ 4 


1.4 x 10~ 5 


Two walls 


1.5 x 10~ 4 


3 x 10" 4 


4 x 10" 8 



TABLE V: Significance of the various fits compared with that 
of the Standard Model — = 0. The smaller a number in the 
Table, the more the model indicated in the corresponding row 
is statistically significant over the Standard Model. Random 
errors are added to the data as indicated in Section [TTT1 



precise solution is not needed for the fit. 



Statistical preference of 
monopole over : 


Keck 


VLT 


Total 


Dipole 


0.60 


0.04 


4 x 10" 5 


One wall 


0.14 


4 x 10" 4 


8 x 10" 5 


Two walls 


0.48 


5 x 10~ 4 


3 x 10~ 7 



TABLE VI: Significance of the various fits with that of a 
monopole template, — = m. Random errors are added to 
the data as indicated in Section [TTT1 



non- vanishing monopole term allows one to study the sig- 
nificance for a spatial variation of a including the possible 
systematic errors discussed in @|. 



VI. SPATIAL DEPENDENCE 

We divide this Section in two parts. In the first part we 
study the significance of the fits for which a exhibits spa- 
tial dependence. In the second part, we further discuss 
the spatial orientation of the different solutions. 



A. Significance of the fits 

For this discussion, we assume that the extra random 
errors included in 0, [HI are a valid estimate of the ac- 
tual error in the data. Here, we pose a well defined ques- 
tion of whether "model 2" (spatial variation), which in- 
cludes "model 1" (no spatial variation) as a sub-case, is 
a significantly better fit to the data. The statistical sig- 
nificance of the improvement is obtained by comparing 

2 _ 2 2 

F = — — — / X2 (where yf and are the y 2 and the 
number of parameters of the i— th model, respectively, 
while n is the total number of observations) of the ac- 
tual data with that of a random set of data. The values 
reported in Table [V] are the probabilities that fitting a 
random set of data with the model indicated in the row 
results in a greater \ 2 improvement over the Standard 
Model than the improvement obtained with the actual 
data. For example, we see that for the VLT data, there 
is a 10% chance that random data would yield a compa- 
rable improvement in \ 2 relative to that of the monopole. 



Table fVTl shows the analogous comparison of the dipole 
template and the wall models against the monopole tem- 
plate. While Ref. [2J discussed a number of uncertainties 
that could affect the value for a obtained for each ab- 
sorber, none of these uncertainties is correlated with the 
spatial position of the absorber. Therefore the uncer- 
tainties studied in Q may result in at most a bias in 
the average value of Aa that one obtains from the data. 
Therefore, comparing the dipole and wall fits against a 



Most of the results presented in these two tables al- 
ready emerged in the discussions done in the previous 
Section. The significance for the monopole over the Stan- 
dard Model (which, as we said, could be interpreted as a 
signal of an overall systematic error in the data) is very 
strong for the Keck absorbers, and less marked for the 
VLT measurements. The Keck data do not show any 
evidence for a dipole modulation being preferred over 
monopole, while a stronger significance emerges for the 
one-wall model. The significance for the dipole is in- 
creased in the VLT data; however, the one-wall solution 
has a far greater significance than that of the dipole. 
The individual VLT and (particularly) Keck measure- 
ments do not support the two-wall model over the one- 
wall model. However, a very strong significance for the 
two- wall model appears once all the data are included. In 
this case, the two-wall model is a significantly better fit 
than the dipole template and the single- wall model, even 
accounting for the fact that it has three more parameters. 



B. Spatial orientation 

Finally, let us further study and compare the spatial 
orientations of the various fits. We confirm that the 
Keck and VLT dipole orientations are consistent with 
each other Q. The orientation of the one- wall fit of the 
Keck data is also in excellent agreement with that of the 
dipole fits as seen in Figure [7] (in the wall fits, the orien- 
tation refers to the line of sight to the point on the wall 
closest to us). This is not the case for the VLT data; 
in this case the dipole and the one-wall orientations are 
only consistent with each other at the ~ 2a level (see 
Figure [8] and Table HII]) . More importantly, the domain 
structures indicated by the Keck and the VLT data are 
incompatible with each other, as they imply an opposite 
sign for A on the other side of their respective walls. The 
inconsistency is also noted by the fact that the best one- 
wall fit of the combined sample is oriented in a different 
direction with respect to the two separate samples (see 
Figure El). 



12 



The separate analyses are indeed rather suggestive of 
the presence of two different domain walls. This is con- 
firmed by the two-wall fits. Neither separate sample pro- 
vides evidence for the two- wall fit: each separate sample 
only covers a limited portion of the sky, with the Keck 
data more concentrated in the northern galactic hemi- 
sphere, and VLT data concentrating on more southern 
galactic coordinates (the different positions of the ab- 
sorbers are clear from comparing Figures [7] and [8] the 
position of the galactic disk in these coordinates can be 
seen in Figure 30 of 0] ) . The evidence for the two domain 
walls only emerges when the two samples are combined. 1 
We remark that, for the total sample of data, the double- 
wall fit is much more significant than either the dipole or 
the single- wall fits. 

In comparing the different fits, one should not be mis- 
guided by the expectation that the one- wall model should 
always be well modeled by a dipole. This is the case only 
if the wall is very close to us, so that most (if not all) 
of the absorbers in the direction of the wall are indeed 
beyond it. Only in this case, the wall model assigns a 
nonvanishing A to those absorbers, and one would find 
that, approximately, the value of a of the absorbers is 
only function of the hemisphere in which they lie, which 
would indeed be well approximated by a dipole; in the 
opposite case, we note that a wall further than all the 
absorbers provides Aa = for all the absorbers, exactly 
like the Standard Model. Therefore, we expect that a 
single model can be well approximated by a dipole only 
if the redshift of the wall is sufficiently small. This is con- 
firmed by the data: only for the Keck sample the best 
dipole fit is a good approximation to the best single-wall 
fit; it is not a coincidence that the value of z wa ii emerging 
from the Keck sample is smaller than the one emerging 
from the VLT and from the combined samples. In this 
regard, the fact that the dipole and the single wall best 
fits of the combined data have a comparable \ 2 appears 
to be just a coincidence. 



VII. SUMMARY 

In this work, we examined models which could account 
for a spatial variation in the electromagnetic fine struc- 
ture constant from Keck and VLT observations. We used 
data derived by 0, H| using the many-multiplet method 
on each absorption spectrum. We fit this set of data 
against one of the dipole templates studied in [2j , against 
the dilatonic domain wall model of (l6j . and against a 
simple generalization of this model containing two walls. 
It is first of all worth stressing that all these models are 



The non uniform distribution of the absorbers on our sky appears 
to be the origin of the non-gaussian x 2 behavior observed in 
Figures [5] and [6] The sky coverage is not sufficiently complete 
and uniform for the central limit theorem to apply. 



very bad fits to the data (see Table A} if one takes only 
the statistical error provided by the many-multiplet pro- 
cedure as the error in the measurements. To compensate 
for this, ref. Q added a common error (in quadrature) to 
27 of the Keck absorbers, and ref. [2[ added a common 
random error to all of the VLT measurements. The value 
of the added errors are chosen such that a monopole (in 
the Keck case) and a dipole (in the VLT case) are good 
fits to the data (we described the procedure in Section 
llllj) . For the Keck data, a case was made for adding such 
an error on 27 objects, since different lines (with respect 
to the other absorbers) are used for deducing the value 
of a in these absorbers. An analogous justification was 
not given for the VLT data. Given the profound con- 
sequences that a variation of a would imply, it would 
certainly be preferable to understand and quantify this 
error before embarking on an analysis of the data. 

Even if this problem precludes us from firmly infer- 
ring a spatial dependence from the data, it is nonetheless 
interesting to analyze the data under the same assump- 
tions made in 0, y], to understand whether the dipole 
template considered in those works is robust. Finding for 
example that some different fit describes the data better 
could help in identifying possible systematic origins of 
the spatial variation. Furthermore, no underlying phys- 
ical model was presented in @, Q for a dipole variation 
of a. The simplest theoretical way to achieve a variation 
of a is to assume that a is controlled by the vacuum ex- 
pectation value of a field; in string theory models this is 
the dilaton field. If the potential of the dilaton admits 
multiple minima, one may imagine a situation in which 
different dilatonic domains exist in the universe, char- 
acterized by different values of a. A variation of a on 
large spatial scales would result if one of these domains 
is at a distance comparable to the typical distances of the 
Keck and VLT absorbers [16( . It is interesting to see how 
this more physically motivated model compares against 
a dipole template. 

As the dipole and the wall fits have the same number 
of parameters (4), an immediate comparison is obtained 
by comparing the resulting \ 2 for each model. The wall 
performs better when the Keck and VLT samples are an- 
alyzed separately. However, the Keck and the VLT best 
dipoles are consistent with each other, while the wall fits 
are not (the Keck wall has A < for the other vacuum, 
while the VLT wall has A > 0). As a consequence, the 
combined data set are fit (marginally) better by a dipole 
than by a single wall; the difference is however of only 
5x 2 — 1.4 over 293 objects in the sample, so that the two 
fits are nearly equivalent. The fact that the two samples 
were best fit by two different walls, and the fact that the 
Keck and the VLT absorbers - with a few exceptions - 
cover different portions of the sky (so that a second wall 
could affect a largely different sample), prompted us to 
investigate the total dataset with two walls. This turns 
out to be the fit with the highest significance, even if it 
has 3 more parameters than either the single wall and 
the dipole fits. With the extra random errors added, the 
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statistical significance of the Standard Model over the 
latter fits is O (10~ 5 ), while the statistical significance 
of the Standard Model fit over the two- wall is O (lO -8 ). 
Even when the extra errors are not included, the two- wall 
model is certainly the best fit of the data. The addition 
of its 7 parameters allow the x 2 to decrease by 100 with 
respect to the Standard Model. When the random errors 
are included, the drop in \ 2 is 47. 

In conclusion, if the procedure adopted in 0, H| for the 
extra errors turns out to be a good measure of the total 
errors, it is clear that the two-wall model gives a viable, 
statistically significant, physical model for the spatial de- 
pendence of the data. The model that we have proposed 
in Section |TT] to account for the multiple walls undergoes 
two stages of symmetry braking as the temperature of 
the universe decreases. Cosmic strings form at the first 
stage (with a tension too smaller to be of any cosmologi- 
cal relevance) , and multiple domain walls stream off each 
string at the second breaking (as discussed in [16| , astro- 
physical and cosmological limits would force the tension 
of the wall to be of ~ O (MeV)). Therefore, the physical 
statement associated with the two wall model is that we 
would happen to live cosmologically close to one of the 



strings formed at the first breaking. On the other hand, 
it is also worth noting that, if the added random errors 
are correct, the Standard Model itself is no longer a ter- 
rible fit to the data. Its value for the total data set 
is 23%. The significances of the dipole and wall models 
over the Standard Model quoted in the previous para- 
graph are purely statistical ones, and refer to an unbi- 
ased model comparison. It is unclear whether we should 
really be unbiased when comparing the Standard Model 
with a template with no underlying model, or even with a 
model that is not supported by any other measurement. 

Note added: During the completion of this work, sim- 
ilar results concerning the position and statistical signif- 
icance of a domain wall in relation to a spatial variation 
in a appeared [29]. 
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